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Cross-Reference to Related Applications 

The present application is related to commonly assigned U.S. patent application Serial No. 

(Attorney Docket No. F1S920030107, titled "EXTENDING THE RANGE OF 

LITHOGRAPHIC SIMULATION INTEGRALS TO INCLUDE INTERMEDIATE DISTANCE 

SCALES"), U.S. patent application Serial No. (Attorney Docket No. FIS920030 1 09, 

titled "IMPROVEMENT OF PERFORMANCE IN SECTOR-BASED OPC ENGINE UTILIZING 
EFFICIENT POLYGON PINNING METHOD AND SIMULTANEOUS COMPUTATION OF 

MULTIPLE SAMPLE POINTS"), U.S. patent application Serial No. . (Attorney 

Docket No. FIS920030 1 1 0, titled "RENESTING FLAREMAP INTO DESIGN FOR EFFICIENT 

LONG-RANGE FLARE CALCULATION"), and U.S. patent application Serial No. _ 

(Attorney Docket No. FIS920030262, titled "PARALLEL COMPUTATION OF MULTIPLE 
POINTS ON ONE OR MULTIPLE CUT LINES"), filed on even date herewith, which are hereby 
incorporated by reference herein in their entirety. 

Background of the Invention 

1. Field of the Invention 

This invention relates generally to the field of optical miciolithography, and more particularly to 
the simulation of a wavefront and the incorporation of a phase map to analytically account for flare 
in optical proximity correction of photolithographic masks. 
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2. Description of Related Art 

The optical microlithography also known as the photolithography process in semiconductor 
fabrication consists of duplicating desired circuit patterns as best as possible onto a semiconductor 
wafer. The desired circuit patterns are typically represented as opaque and transparent regions on 
a template commonly called a photomask. In optical microlithography, patterns on the photomask 
template are projected onto photoresist-coated wafers by way of optical imaging through an 
exposure system. 



Aerial image simulators, which compute the images generated by optical projection systems, have 
1 0 proven to be a valuable tool to analyze and improve the state-of-the-art in optical lithography for 
integrated circuit fabrication. These simulations have found appHcation in advanced mask design, 
such as phase shifting mask (PSM) design, optical proximity correction (OPC), and in the design 
of projection optics. Modeling aerial images is a crucial component of semiconductor 
manufacturing. Since present lithographic tools employ partially coherent illumination, such 
1 5 modeling is computationally intensive for all but elementary patterns. The aerial image produced 
by the mask, i.e., the light intensity in an optical projection system's image plane, is a critically 
important quantity in microlithography for governing how well a developed photoresist structure 
replicates a mask design. 

20 Optical proximity correction simulation kernels associated with lithographic processes for 
semiconductor chip manufacturing currently do not take into account the higher order aberrations, 
whose resuhs may be realized in long-range effects, although not as prominent in the close-range 
of 1 to 2 micrometer. As the state of the art moves towards smaller wavelengths of light, such as 
193nm and 157nm and extreme ultraviolet (EUV) 13nm; and with device dimensions becoming 

25 considerably smaller in ratio to the wavelength of light that is used to print them on the wafer, the 
long-range effects, such as flare, become significant, making it imperative that higher order 
aberrations be considered. 

In the prior art, data analyses have been hampered by the lack of a flare-capable simulator for 
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handling large areas of a mask. Flare is generally defined as unwanted light in a lithographic process 
located in places where it should otherwise be dark. The cause of this is threefold. First, waveftont 
roughness or high order wavefront aberrations, which encompass up to approximately 10'° 
Zemikes, cause flare. Optics polishing limitations, contamination, and index variations, such as 
frozen turbulence inside fused silica, all contribute to wavefront roughness. Second, ghost images 
or stray reflections, on the order of 1% for a 99% antireflective coating, will also contribute to flare. 
Last, the scattering from the walls of projection optics may cause flare. 

If the flare were constant, a dose shift would compensate for its effects completely. However, it is 
not constant, and a 1% flare variation may result in an unacceptable Across Chip Line-width 
Variation (ACLV). ACLV is a key parameter used to describe the fidelity of the printing process. 
In a design line-width of 45mn, ACLV of 6mn represents about +/- 1 5 % variation. The acceptable 
variation is nominally +/-10% or less. This large variation can considerably diminish the circuit 
performance and in some cases may cause a catastrophic failure. Thus, it is necessary to determine 
and compensate for the flare's contribution to the ACLV. 

MathematicaUy, flare represents the sum of all halo contributions from the bright regions of the 
mask; the light surrounding theoptical system point spread fimction(PSF) that is causedbyscatter 

from within the optical system. Thus, flare scales with the bright area, hi the short range, which is 
estimated to be on the order of R,^, ~ 5Xy(2NA) out to 2.5 micrometers (whereNA is the numerical 
aperture of the optical system and X represents the wave length of the light) partial coherent effects 
are realized. Flare contributions add up incoherently in the medium range (2.5 to 5 micrometers) 
and m the long range (5 micrometers up to ~ 10 mm). 

Other long-range effects include non-optical effects such as etch, macro loading effects, and 
chemical flare. 

Flare increases with the areaofbrightness. Thus, bright field masks, where thebackgroundisclear 
or bright, are susceptible to flare in general. Double exposure, where two exposures with two 
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different masks are used to print a single set of shapes, may also cause more of a flare problem. 
Light effectively leaks from the shield during the second exposure, adding new flare to that already 
produced in the primary exposure. Masks with varying pattern densities are more susceptible to 
flare variation. Under certain physical conditions, flare scales as the inverse square of the 
wavelength, 1 /X^ and canbecome more problematic in smaller wavelengths of light including EUV. 

Experiments have demonstrated complex flare effects at multiple scales. For example, kernel 
anisotropy, field variations, and chemical flare have all been shown. Until now, however, data 
analysis has been hampered by the absence of a flare-capable simulator. 

In U.S. Patent No. 6,263,299 issued to Aleshin, et al., on July 17, 2001, entitled, "GEOMETRIC 
AERIAL IMAGE SIMULATION," an aerial image produced by a mask having transmissive 
portions was simulated by dividing the transmissive portions of the mask into primitive elements, 
obtaining a response for each of the primitive elements, and then simulating the aerial image by 
combining the responses over all of the primitive elements. However, this analytical approach does 
not teach or suggest a method to account for the higher order aberrations of flare within an optical 
proximity correction simulation kernel. 

In prior art, the shape of the wavefront is represented by a series whose terms are orthogonal. The 
most commonly used terms are the Zemike polynomials. To date, higher order aberration effects 
have not been analytically expressed as polynomials in the art, in contrast to the expression of 
Zemike polynomials for lower order terms. In the current state of the art Optical Proximity 
Correction (OPC) tools, the number of Zemike polynomials used is limited up to on the order of 
37. Consequently, there exists a need in the art to provide a method of accommodating the higher 
order aberration effects in subsequent model calibrations. 

Accommodating higher order aberrations will ultimately allow for: accuracy in the computation of 
calibrated optical and resist models; fidelity in the wafer shapes to the "as intended" shapes by 
achieving better correction of mask shapes; accuracy in simulation of wafer shapes for a better 
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understanding and evaluation of correction methodologies; increases in yield in chip manufacturing 
due to the better accuracy; and reduced cost in fabrication. 



Summary of the Invention 

Bearing in mind the problems and deficiencies of the prior art, it is therefore an object of the present 
invention to provide a method to account for analytical wavefront construction over the mask at all 
distances. 

It is another object of the present invention to provide a method for creating a phase map from the 
wavefront. 

It is a further object of the present invention to provide wavefront-based kemels for implementation 
in the short, long, and medium range optical effects. 

It is another object of the invention to provide a method of accommodating the higher order 
aberration effects utilizing existing manufacturer's data for wavefront construction. 

It is yet another object of the present invention to provide a method of accommodating the higher 
order aberration effects utilizing simulated wavefront data. 

Still other objects and advantages of the invention will in part be obvious and will in part be 
apparent from the specification. 

The above and other objects, which will be apparent to those skilled in art, are achieved in the 
present invention, which is directed to a method of performing model-based optical proximity 
correction comprising: embedding wavefront information on a first two-dimensional complex array 
having a plurality of array elements and an assigned diameter; generating a phase map from the 
wavefront information; computing a point spread function from the phase map; and performing 
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optical proximity correction calculations using the point spread function. The step of embedding 
wavefront information comprises embedding simulated wavefront information from randomly 
generated data, or embedding empirically derived wavefront data. 

The method further comprises: forming the first two-dimensional complex array comprising a 
plurality of complex numbers having amplitude and phase; obtaining a complex conjugate array 
corresponding to the first two-dimensional complex array; arranging the first two-dimensional 
complex array into a symmetric complex array using the complex conjugate array; performing an 
analytic filtering fimction on the symmetric complex array; applying a power law fiinction to the 
symmetric complex array; computing the wavefront by performmg a Fourier Transform on the 
symmetric complex array, resulting in a transformed array; obtaining a circular core of the 
transformed array of a diameter equal to the first arrays assigned diameter, and converting to zero 
coordinates outside the diameter; and normalizing amplitudes of the transformed array. 

The step of fonning the first two-dimensional complex array comprises putting the array elements 
in phasor form with a phase within a range of 2p. The phasor may have an amplitude in the range 
of 0 to a predetermined small fraction less than 1 . The first two-dimensional complex array size can 
be represented by the log in base two of the first two-dimensional complex array's assigned 
diameter. The complex numbers are centered within the first two-dimensional complex array, 
equatingeach amplitudeofthefirst two-dimensional complex array withacorresponding amplitude 

of the complex conjugate array, and equating each phase of the first two-dimensional complex array 
coordinate pair with the inverse of a corresponding phase of the complex conjugate array. 

The step of applying a power law fiinction on the symmetric complex array comprises multiplying 
the amplitude of each of the complex number by the sum of the squares of each of the coordinate 

pair, the sum raised to a power. 

The step of normalizing the transformed array's amplitudes comprises: subtracting an average value 
of the wavefront from each array element; and replacing each of the elements by an analytical 
expression, which is a fiinction of both a square root of intrinsic flare and a radius of the wavefront. 
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The method further includes shuffling the elements such that the element in the first array's center 
is shifted to the first array's bottom-left comer. This is accomplished by: ignoring a first row and 
first column of the first array; dividing the first array's remaining portion into four quadrants; 
exchanging a first quadrant with a second quadrant, and exchanging a third quadrant with a fourth 
quadrant; exchanging the first and the fourth quadrants; exchanging the second and the third 
quadrants; converting the elements into real and imaginary form; unshuffling the elements such that 
the first array is again centered following the Fourier Transform; and converting the first array into 
phasor form. 

The method may also include: exponentiating empirically derived wavefi-ont information resulting 
in an array of complex numbers in phasor form; embedding the exponentiated wavefront in a larger 
array having a center such that if a Unear dimension of the wavefiront array is given by log diameter 
p, then a corresponding linear dimension of the larger array is given by log diameter q, where q is 
at least p+3; shuffling the larger array to move the center to a bottom-left quadrant; performing an 
inverse Fast Fourier Transform on the larger array to obtain a transformed array; unshuffling the 
transformed array to move the complex numbers fi-om the bottom-left quadrant back to the center; 
normalizing the magnitude of the unshuffled array; choosing amplitude values of the normalized 
array and creating a phase map; using the phase map to compute a point spread function array; and 
scaling the point spread fianction array fi-om pixel dimension to real dimension. 

In a second aspect, the present invention is directed to a method of performing model based optical 
proximity correction on a lithographic mask pattern incorporating phase maps comprising: 
incorporating a point spread function array in real dimension within a set of convolution kernels; 
and computing an aerial image with aberrations using the set of convolution kernels. 

In a third aspect, the present invention is directed to a program storage device readable by a 
machine, tangibly embodying a program of instructions executable by the machine to perform 
method steps for performing model-based optical proximity correction, the method steps 
comprising: embedding wavefi-ont information on a first two-dimensional complex array having 
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a plurality of array elements and an assigned diameter; generating a phase map from the wavefront 
information; computing a point spread function from the phase map; and performing optical 
proximity correction calculations using the point spread function. 

The method further comprises: forming the first two-dimensional complex array comprising a 
plurality of complex numbers having amplitude and phase; obtaining a complex conjugate array 
corresponding to the first two-dimensional complex array; arranging the first two-dimensional 
complex array into a symmetric complex array using the complex conjugate array; performing an 
analytic filtering fiinction on the symmetric complex array; applying a power law fiinction to the 
symmetric complex array; computing the wavefront by performing a Fourier Transform on the 
symmetric complex array, resulting in a transformed array; obtaining a circular core of the 
transformed array of a diameter equal to the first array's assigned diameter, and converting to zero 
coordinates outside the diameter; and normalizing amplitudes of the transformed array. 

Brief Description of the Drawings 

The features of the invention believed to be novel and the elements characteristic of the invention 
are set forth with particularity in the appended claims. The figures are for illustration purposes only 
and are not drawn to scale. The invention itself, however, both as to organization and method of 
operation, may best be understood by reference to the detailed description which follows taken in 
conjunction with the accompanying drawings in which: 

Fig. 1 depicts a wavefront roughness relationship to flare. 

Fig, 2 is a graph of the power spectral density of a wavefront. 

Fig. 3 depicts the segments of the incoherent integration scheme as a fiinction of shape density on 
the mask. 
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Fig. 4 depicts the overall architecture of an algorithm for implementing the computation of an aerial 
image at a point that incorporates the wavefront. 

Figs. 5A-5C compare simulations of the wavefront derived from the power law calculations with 
real wavefront. 

Figs. 6A and 6B depict the point spread function with flare for the gamma values 1.7 and 2.3, 
respectively, as a function of range R. 

Fig. 7 is a flow chart of the process steps that analytically transform a random array A into real 
wavefront data. 

Fig. 8 depicts a resultant simulated wavefront as created by the methodology of the present 
invention. 

Fig. 9 depicts the process flow for the method of embedding measured or empirically derived data 
to obtain a real wavefront array. 

Fig. 1 0 depicts a picture of a wavefront created from a file of measured or empirically derived data. 
Fig. 1 1 depicts an example system capable of implementing the present invention. 
Description of the Preferred Embodiment(s) 

In describing the preferred embodiment of the present invention, reference will be made herein to 
Figs. 1-11 of the drawings in which like numerals refer to like features of the invention. 

The present invention provides a method to compute a phase map that can be used to take into 
accountthehigher order aberrations within an optical proximity correctionsimulationker^ 
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results in a more accurate computation of OPC when long-range effects become prominent. 

Inverse Fourier transforms are used to create a phase map that accounts for higher order aberrations 
normally ignored in traditional optical proximity correction (OPC) simulation kernels. Two 
methodologies are proposed to account for the higher order aberrations. The first method utilizes 
simulated wavefront infomiation from randomly generated data. The second method uses measured 
or empirically derived data from optical tools. 

In both methods, flare is accounted for. Fig. 1 depicts a wavefront roughness relationship to flare. 
A single frequency wavefront 12 is shown with zero and first order transmissions 14 on a wafer. 
The intensity pattern 16 on the wafer is given in square amplitudes and shown as a fimction ofthe 
number of cycles per pupil, the wavelength, and the numerical aperture. Real wavefronts have a 
continuous disfribution of frequencies. The power spectral density (PSD) ofthe transmitted light 
can be calculated from the square ofthe absolute value ofthe Fast Fourier Transform (FFT) ofthe 
wavefront using the following relationships: 
psd ~ |FFT of wavefrontf ~ Amp^ 
and 

Amp' ~ Power Law Function ~ l/(#cycles/pupiiy ~ 1/x' - Flare Halo 
where, 

#cycles/pupil ~ (Zemike #)"' 

Polishing and index fluctuations yield y values of approximately 2. However, the flare halo rapidly 
exceeds diffraction, creating a diffraction-limited image. 

The flare is a fimction ofthe intrinsic flare and the integral of the power specfral density and mask 
transmission. When the mask transmission is assumed to be 1, the intrinsic flare value is (2na)^ 
where a is the real residual rms wavefront. Fig. 2 is a graph of the power spectral density of a 
wavefront. This figure also denotes the limitations ofthe usage of Zemike polynomials and the 
diffraction limited Point Spread Function (PSF) in simulating flare. The logarithm ofthe power 
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spectral density (PSD) 20 and also the log of the flare halo are graphed on the ordinate axis as a 
function of the number of cycles per pupil 24 in case of (PSD), and as the log of the distance scale 
in case of the log of Halo. The limitation posed by using only on the order of 37 Zemike 
polynomials is depicted at the x-axis (22). Curve A depicts the area for non-optimal interferometers, 
which fundamentally filter out high frequency data. Curve B depicts the area for diffraction limited 
point spread function. This depicts the limitation of the current state of the art OPC tools that use 
only diffraction limited PSF. Curve C depicts a real point spread function, showing a slope of -y. 
Line 26 represents the pixel cut-off for diffraction limited PSF. At this cut-off point, there is enough 
analytical range to determine y. 

Flare at any given point is the sum of the contributions from all other points. Under the integration 
scheme of the present invention, the summation grid may be coarse at certain distant parts of the 
mask. Fig. 3 depicts the segments of the incoherent integration scheme 30 as a function of shape 
density on the lithographic mask. This can also be perceived as the intensity density on the wafer. 
Polygon integration 32 is performed on real estate on the order of 10 square micrometers on the 
chip, while partial coherent effects 34 are determined at approximately 3 square micrometers. Note 
that the range is depicted in square for the illustrative purpose only. The ranges may as well be 
described by circles, or any other closed, and convex geometric shapes. Also the distance scales are 
described for illustrative purpose only. They may vary within the same range on order, depending 
on the lithographic process and technology. 

In the intermediate range of approximately 1 .5/^m to S^m (that is within the depicted square of 10 
Aim), the algorithm of the present invention speeds calculations by integrating directly over the 
mask polygons. This avoids summation over every pixel when the required pixel-grid is very fine. 
Additionally, incoherent convolution using an averaged point spread fimction is accurate for large 
polygons at an appreciable distance from the image point. However, edges close to each other or 
to the image point interact coherently. The algorithm of the present invention allows these 
limitations to be partly relaxed. A coherent integration is used inside the short-range zone of less 
than approximately 1 .S/zm (that is within the depicted square of 3^ml and an exact, non-average 
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point spread function is utilized in the intermediate range for incoherent light, or alternatively, the 
average point spread function and analytic solution may be used. Importantly, in order to use and 
exploit these features, the wavefront/point spread function must be known or the user must employ 
simulated wavefronts generated to obey user-input statistics, 

5 

Fig. 4 depicts the overall architecture of an algorithm for implementing the computation of an aerial 
image at a point. First, model options 42 are input by the user. This data is used for short range, 
partially coherent and exact geometry calculations 44 to obtain the short-range effects of the aerial 
image 46. It may also be used to determine wavefront data 48. Empirically derived wavefront data 

10 or simulated wavefront data is employed and a point spread function computed 50 for the partially 
coherent, exact geometry, mid-range calculations 52. If no wavefront data is used, the effects of 
incoherent light may be calculated 54. In either case, the mid-range aerial image is then computed 
56 and added to the short-range calculation result 58. Last, the model inputs are used to compute 
the density map for the entire chip 60. A power law kemel is then convolved with fuzzy geometry 

15 62 to calculate the long-range aerial image 64, which is then added to the short and mid-range 
calculation results 58. 

The wavefront technique of the present invention is represented in the overall algorithm of Fig. 4 
by the Simulate/read Wavefront and Compute Point Spread Function task 50, and the calculating 
20 of the SOCS, partially coherent, exact geometry mid-range effects 52. Figs. 5A-5C contrast images 
of the wavefront derived from the power law calculations 5 A, 5B to a measured wavefront 5C. The 
power law estimates are shown for gamma=1.7. Fig. 5A, and gamma = 2.3, Fig. 5B. Note that a 
smaller gamma value denotes more high-frequency terms for the wavefront as shown in Figs, 5 A 
and5B. 



25 



Figs, 6A and 6B depict the point spread function with flare for the gamma values 1.7 and 2.3, 
respectively, as a fimction of range R. Curves A and A represent a 10% flare, while curves B and 
B' represent a 0% flare. Note again that a smaller value of gamma shows a wider spread of the PSF 
denoting a wider effect for the flare. This is shown in Figs. 6 A and 6B. 
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The simulated wavefront method concentrates only on the high frequency components of the 
aberrations while ignoring the low frequency components by using an apodizing function such as 
a super-Gaussian function. This is because the lower order aberration terms are usually included 
as the Zemike polynomials in the OPC computations. Essentially, the wavefront is filtered by a 
Super Gaussian function to attenuate the low frequency contributions to the wavefront, while 
allowing the high frequency contributions to remain intact. Note that the any other attenuating 
function that effectively filters out the lower order frequency terms can be used instead of the Super 
Gaussian function. An inverse power law fimction is applied on the frequency components so that 
the fiinction emulates the current experimental results of the long-range effects of the optics. The 
wavefront is generated from this information by applying a Fast Fourier Transform. The resulting 
wavefront frinction is then made real by constructing a symmetric complex array A that satisfies 
the equation A(x,y) = A*(-x,-y). 

A phase map is created from either a simulated wavefront or an empirically derived wavefront by 
analytically embedding the wavefront on an array having a significant number of guard band zeros, 
and performing an inverse Fourier transform on the resuhant array. Following the inverse Fourier 
transform, the pixel dimension is then converted to the nominal dimension for implementation with 
OPC kernels. 

The method to generate the phase map is the same whether the wavefront is obtained from 
simulated data or empirically derived information. Simply described, the methodology comprises 
either the generation of a simulated wavefront or the embedding of measured wavefront data, and 
the generation of a point spread fiinction of the wavefront, while retaining the higher order 
contributions. 

Methodology to Generate a Simulated Wavefront 

The methodology for creating a simulated wavefront begins with a randomly generated array and 
ultimately ends with the output of a real wavefront. First, a randomly generated two-dimensional 
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array A of complex numbers is constructed. The elements of array A may be given in phasor form, 
such that phase and amplitude information is retained. This may also be accomplished by building 
a complex array of real and imaginary coordinates. The phase is made to vary uniformly in a range 
of 27C, preferably from -n to +71, and the amplitude made to vary from 0 to a user specified number, 
which may, for example, be on the order of 0.01 . 

The elements of array A are arranged such that for any given (x,y) coordinate, the value of array 
element A(x,y) is equal to the conjugate of the value at the inverse of the coordinates. This may be 
expressed as the following relation: 

A(x,y) = A*(-x,-y), 
Where A* is the conjugate of A. 

The conjugate array must satisfy an equality relationship for amplitude and an inverse equality 
relationship for phase: 

amplitude(A) = amplitude(A*); and 
phase(A) = - phase(A*). 

The (x,y) coordinates are assigned so that they are centered within array A. To achieve this, the 
numbering of the rows of array A are varied from - y_max to +y_max and numbering of the 
columns of array A are varied from -x_max to +x_max, leaving the origin at the array center. This 
numbering is the preferred scheme; however, any other numbering scheme having a similar effect 
may be used. 

To increase the efficiency of the FFT software used, we represent the array size as a number that 
is something very close to a power of 2. However, to ensure the symmetry of A so that it satisfies 
the equation A(x,y) = A*(-x,-y), the number of rows and columns of A needs to be an odd number. 
Therefore, the number of rows and columns of A is expressed as 2^+1, where p is referred as the 
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log-diameter of A. This is because p is equivalent log (#of rows of A -1). Also, the size of A may 
then be determined by the log of the array's diameter, which is an integer. For example, for an initial 
input size of 7, the size of A is given as 129 x 129, where the length 129 is determined from the 
relation 2^+1. 

Once the randomly generated two-dimensional array A of complex numbers is constructed, the 
lower order Zemike terms are filtered out. This is performed so that the terms are not considered 
twice in the same calculation. The filtering requires the amplitude of each element of A to be 
multiplied by a Super Gaussian function or other appropriate filtering fimction. The preferred 
fimction is a super Gaussian fiinction, sg(z). The super Gaussian function is determined from the 
following relationships: 

f=X./(R_min*NA), 

where R_min is the radius of the circular region for the lower order 
terms; 

NA is the numerical aperture; and 
X is the wavelength of light; 

and 

p = ln(l - ln(9)/ln(10))/(2 ln(l - f)) 

The super Gaussian function sg(z) is then defined from these relationships as: 

sg(z) = exp{-ln(2)||(zf)P||} 

In order to achieve the necessary filtering, the amplitude of each element of array A(x,y) is replaced 
by the super Gaussian function acting on the amplitude, such that the resultant amplitude equals 

sg(v/(xV)). 

A power law is then applied to array A by multiplying the resultant amplitude of each element of 
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array A by {{x^-\- ff^^^}, where y is a user supplied function preferably of a value varying from 1 
to 3. 



The wavefront is then computed by applying a Fourier transformation to filtered array A. A 
shuffling mechanism may be applied to the array to enable the array to be used in a commercially 
available Fast Fourier Transformation tool. This shuffling mechanism transforms the array suitably 
so that the center of the array moves to the lower left hand comer and all other elements move 
appropriately. However, this may not be necessary if the Fast Fourier Transform (FFT) method is 
equipped to handle an array as described above. The shuffling mechanism is described below. First, 
the elements of array A are shuffled so that the origin of array A (which was at the center of A) is 
relocated at the bottom-left comer of array A. This is done by the following steps: 

i) ignore the first row and the first column of A; 

ii) divide the elements of A in four quadrants; 

iii) exchange the 1 st and the 2nd quadrants; 

iv) exchange the 3rd and the 4th quadrants; 

v) exchange the 1 st and the 4th quadrants; and 

vi) exchange the 2nd and the 3rd quadrants. 

Note that step (i) causes the number of rows and columns in the array that is passed to the FFT to 
be a power of 2, if the number of rows and columns of A is initially set at 2p+1. The elements of 
the shuffled array A are then converted from the phasor form to real/imaginary form. This is done 
because the current commercially available FFT tools handle data in the real/imaginary form. This 
may not be necessary if the FFT tool can handle the data in the phasor form, A two-dimensional fast 
Fourier transform is then applied to the array. The array is then unshuffled so that it is again 
properly centered. Next, it is converted back to phasor form. 

After the Fast Fourier Transform (FFT), the phase of the elements of array A will be zero or 2p. 
This is ensured by the initial assumption that the array coordinates equal the complex conjugate of 
the inverse coordinates, A(x,y) = A*(-x,-y). 



FIS920030108US1 



-16- 



A circular core of array A is then taken. The diameter of array A is same as the dimension of the 
original array; however, the rest of the elements are converted to zero. The amplitudes are then 
normalized in the following manner: 

i) the average value of the amplitude is subtracted from each element; 
5 ii) amplitude of each element A(x,y) is replaced by the expression: 

(^(Intrinsic Flare)/47i) / (v/(S(A(x,yy)/ n . r), 

where "r" represents the radius of the wavefront. 

The resultant wavefront is real because the phase of all the elements has a value of zero or In, Fig. 

10 7 is a flow chart of the process steps 80 that analytically transform a random array A 82 to real 
wavefront data 94. As previously discussed, the random array A is made to satisfy the condition 
A(x,y) = A*(-x,-y), 82. The array core is then apodized to remove the lower order Zemikes 84. A 
power law function is applied 86. The array is then shuffled, subject to a fast Fourier transform, and 
unshuffled 88. The core is circularized 90 and normalized 92, resulting in a wavefront represented 

1 5 by real data. Fig. 8 depicts a simulated wavefront as generated by the above-described methodology. 

Method for Embedding Measured Wavefront Data 

The second methodology proposed to account for the higher order aberrations uses measured or 
20 empirically derived data from optical tools. The wavefront data is usually available from tool 
manufacturers in the form of raw data. This data needs to be properly embedded and centered 
before it can be used to compute the point spread frmction. The data provides the following 
information: 

i) the raw data; 

25 ii) a description of values that can be analytically ignored; and 

iii) the nimiber of rows and number of columns, which need not be same. 
The following steps are employed to embed the data: 

i) reading the data in row major order by substituting a zero for each ignored piece of 

data; 
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ii) centering the data, which is done to avoid unintentional tilt in the data in the 
following manner, as follows: 

a) locating the best circle that encloses the data; 

b) embedding and centering the circle within a square array, whose dimension is 
given as the log diameter, p, such that the number of rows and columns of the 
array is given as 2P^'. 

Fig. 9 depicts the process flow 100 for the method of embedding measured or empirically derived 
data to obtain a real wavefront array. As shown in the figure, the measured data is inputted 102, 
circularized with the best fit circle 104, embedded and centered 106, fitted in a square array with 
the best fit log-diameter 108, and converted to real wavefiront data 1 10. 

A picture of the wavefiront created firom a file of measured data is shown in Fig. 10. 

Methodology to Create a Point Spread Function firom Real- Valued Wavefront Data 

Whether the wavefi-ont is generated from simulated or empirically derived data, a point spread 
fiinction is created that will be used to compute the aerial image that accounts for higher order 
aberrations. First, the real-valued wavefront W is exponentiated as e'^. Note that the wavefront 
created so far is a real array with phase values being zero. Therefore exponentiating this real array 
has the resulting effect of an array of complex numbers in the phasor form (phase and amplitude), 
whose phase values correspond to the amplitude values of the wavefront (this follows from the 
definition of exponentiation), and the amplitude values of all the elements are unity. The wavefront 
is then embedded in a larger array. As an example, if the dimension of the wavefront array is given 
by a log-diameter p, then the dimension of the larger array is given by a log-diameter q, where q is 
preferably at least p+3. Note that for an array, if the dimension is given in terms of a log-diameter 
p, then the number of rows and columns of the array is given as 2^+1 . 

Next, as before, we prepare the array for an inverse FFT by using a shuffle. This step may not be 
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needed if the FFT tool can handle the data in the given fomat. The array is shuffled to move the 
center of the array to the bottom left. An inverse fast Fourier transformation is performed on the 
shuffled array. The array is then unshuffled to move the bottom-left portion back to the center. The 
array is scaled from the pixel dimension to the real dimension. 

5 

The above-described method computes a phase map that can be used to take into account the higher 
order aberrations within an optical proximity correction simulation kernel. The method may use 
either simulated wavefront information or empirically derived information. The wavefront 
information is modeled within an array, whose elements are reordered, given a complex conjugate 
1 0 correlation, operated upon by a filtering function, and subjected to a fast Fourier transform and later 
an inverse Fourier transform. The resultant wavefront is calculated to demonstrate the affects of the 
higher order aberrations. 

The method may be implemented through software in a program storage device readable by a 
1 5 machine, tangibly embodying a program of instructions executable by the machine to perform each 
of the method steps. Fig. 1 1 depicts an example system capable of implementing the present 
invention. CPU 1001, I/O device 1002 and conmiunications device 1003 connect via a bus 1005 
to a central computer system 1004 having memory banks 1006 to store the program instructions. 
Other hardware/software combinations are feasible, and can be configured to provide the necessary 
20 computing power to perform the present method. 

While the present invention has been particularly described, in conjunction with a specific preferred 
embodiment, it is evident that many alternatives, modifications and variations will be apparent to 
those skilled in the art in light of the foregoing description. It is therefore contemplated that the 
25 appended claims will embrace any such alternatives, modifications and variations as falling within 
the true scope and spirit of the present invention. 

Thus, having described the invention, what is claimed is: 
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